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Abstract 

Self-diffusion  coefficients  of  several  different  solvents  in  polystyrene,  as 
well  as,  methyl  methacrylate  in  poly(methyl  methacrylate)  have  been  determined 
at  polymer  concentrations  from  0  to  50  wt%  at  25°C.  The  solvents  used  for  the 
polystyrene  solutions  were  toluene,  ethylbenzene,  cumene,  t-butyl  acetate, 
chloroform,  and  methyl  ethyl  ketone.  The  diffusion  data  was  used  to  evaluate 
several  theories  for  the  concentration  dependence  of  the  solvent  diffusion 
coefficients,  and  in  most  cases  the  reduced  diffusion  coefficients,  of  solvents  in 
pol)niier-solvent  systems.  In  addition  to  the  experimental  data  presented  here, 
literature  data  for  several  systems  were  also  used  to  test  the  models.  The  theories 
evaluated  were  based  on  widely  different  concepts  including  electrical 
conductivity,  tortuosity,  kinetic  theory,  and  free- volume.  The  results  indicate  that 
the  kinetic  theory  and  free-volume  approaches  fit  the  data  most  consistently,  with 
the  kinetic  theory  approach  pieferred  at  lower  polymer  concentrations  and  the 
free-volume  approach  favored  at  higher  polymer  concentrations. 
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Introductimi 

DifTusion  of  small  molecules  in  polymer  solutions  is  often  important  in  a 
variety  of  processes  including:  polymer  processing, polymerization  kinetics, 
drug  delivery,®  and  drying  of  coatings.^  It  would  be  useful  in  many  of  these  areas 
if  the  diffusion  coefficient  of  a  particular  molecule  at  a  particular  concentration 
and  temperature  could  be  determined  easily  without  having  to  measure  the  value 
experimentally.  We  have  found,  to  a  first  approximation,  the  concentration 
dependence  in  the  0  -  50%  polymer  range  can  be  predicted  from  knowledge  of  the 
diffusion  coefficient  of  the  pure  solvent®*^*^®  for  single-solvent  polymer  systems. 
However,  for  more  detailed  knowledge  of  these  dependencies,  or  use  at  higher 
concentrations,  experimental  measurements  appear  to  be  required.^®  To  help 
describe  solvent  diffusion,  models  based  on  electric  conductivity  by  MaxwelH^ 
and  Fricke,^^  and  on  tortuosity  in  a  simple  cubic  lattice  model  by  Mackie  and 
Meares^®  have  been  used^®’^^’^^  because  they  offer  a  simple  dependence  of  the 
reduced  diffusion  coefficient  on  the  polymer  volume  fraction.  While  such  a 
simple  dependence  is  sought,  these  approaches  do  not  account  for  subtle 
differences  observed  in  experimental  data.  In  this  paper,  we  will  explore  two 
models  which  much  more  accurately  describe  the  diffusion  data. 

The  first  model  is  based  on  the  modified  Enskog  theory,^®  which  uses 
kinetic  theory  and  accounts  for  momentum  transfer  between  rigid  spheres.  This 
theory  gives  a  rather  involved  dependence  of  the  solvent  diffusion  coefficient  on 
concentration.  The  second  theory  can  be  used  to  predict  the  dependence  of  the 
reduced  diffusion  coefficient  on  polymer  weight  fraction  from  Vrentas-Duda  free- 
volume  theory. This  results  in  a  simple  expression  with  one  fitable 
parameter  for  systems  in  which  a  variety  of  parameters  are  known.  We  will 
examine  the  effectiveness  of  these  two  theories  for  describing  the  concentration 
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dependence  of  solvent  diffusions  coefficient  by  fitting  the  proposed  equations  to  our 
own  experimental  data,  as  well  as,  experimental  data  taken  from  the  literature. 
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Experimental 

Pulsed'Gradient  Spin-Echo  NMR.  The  diffusion  measurements  were 
performed  using  the  PGSE-NMR  technique.^>^®’^®  In  this  technique,  a  spin-echo 
experiment  is  performed  and  in  between  the  radio  frequency  (RF)  pulses, 
magnetic  field  gradient  (MFG)  pulses  are  applies.  The  MFG  pulses  result  in  an 
attenuation  of  the  echo  amplitude.  The  resulting  echo  amplitude  is  given  by  the 
relationship 


A  =  Ao  exp(-t/T2)  expf-y^G^DP)  ( 1) 

where  A  is  the  echo  amplitude,  x  is  the  time  between  the  RF  pulses,  T2  is  the  spin- 
spin  relaxation  time,  y  is  the  magnetogyric  ratio  which  is  constant  for  a  given 
type  of  nucleus,  G  is  the  strength  of  the  magnetic  field  gradient,  D  is  the  self¬ 
diffusion  coefficient,  and  p  =  6^(A-6/3),  5  being  the  width  of  the  gradient  pulse  and 
A  the  time  from  the  start  of  first  gradient  pulse  to  the  second.  If  an  array  of 
experiments  with  varying  values  of  5  are  done,  the  self-diffusion  coefficient  of 
each  component  can  be  determined  as  long  as  a  unique  resonance,  for  that 
component,  can  be  resolved.  This  is  accomplished  by  plotting  the  natural 
logarithm  of  the  amplitude  of  a  resonance  for  the  compound  in  question  versus  p. 
The  slope  of  the  resulting  line  is  proportional  to  its  diffusion  coefficient.  Since  the 
slope  of  this  line  is  also  proportional  to  y  the  most  convenient  nuclei  (largest  y)  for 
diffusion  studies  are  protons  which  were  observed.  Experiments  were  performed 
on  a  Varian  VXR-200  which  we  have  modified  for  diffusion  studies.^®’^^*^^  Typical 
gradients  used  were  0.16  T/m,  with  pulse  widths  (5)  up  to  5  ms  and  echo  times 
(2x)  of  200  ms.  Often  the  pol3mier  and  solvent  resonances  did  not  overlap,  so  that 
the  solvent  self-diffusion  coefficients  were  easily  obtained.  In  the  case  of  overlap. 
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echo  times  200  ms  were  sufficient  to  remove  the  polymer  signals  from  the  spin- 
echo  spectra  due  to  the  relatively  short  T2  relaxation  times  of  the  polymers. 

Materials.  The  PS  (Aldrich  Chemical  Company,  Milwaukee,  WI)  with  a 
molecular  weight  of  280,000  g/mol  was  used.  The  solvents  used  were  t-butyl 
acetate(Aldrich);  chloroform,  methyl  ethyl  ketone,  and  toluene  (Fisher  Scientific, 
Fair  Lawn,  NJ);  and  ethylbenzene  and  cumene  (Eastman  Kodak  Company, 
Rochester,  NY).  Toluene,  ethylbenzene,  and  cumene  are  similar  chemically; 
while  chloroform,  methyl  ethyl  ketone  (MEK)  and  t-butyl  acetate  vary  in  size,  as 
well  as,  polarity.  To  ensure  that  the  observed  behavior  was  not  limited  to  PS 
systems,  the  diffusion  of  methyl  methacrylate  (MMA)  in  poly(methyl 
methacrylate)  (PMMA)  was  also  studied  (both  from  Aldrich).  The  molecular 
weight  of  the  PMMA  was  12,500  and  polymer  concentrations  between  0  and  50% 
were  studied  for  these  systems.  All  materials  were  used  as  received.  All  of  the 
experimental  results  reported  for  these  systems  are  the  average  of  multiple  (two 
to  five)  runs. 

Literature  Data.  To  broaden  the  scope  of  this  study  data  was  also  taken 
from  the  literature  and  fit  with  the  two  theories.  The  literature  data  was  for 
methyl  ethyl  ketone  in  PMMA  (M.W.  960,000),^^  methyl  ethyl  ketone  in  poly(n- 
butyl  methacrylate)  (PBMA)  (M.W.  450,000),^^  toluene  in  poly( vinyl  acetate)  (PVA) 
(M.W.  500,000),^®  tetrahydrofuran  in  PS  (M.W.  540,000),  and  tetrahydrofuran  in 
polyisoprene  (PIP)  (M.W.  253,000).^®  The  last  three  systems  were  studied  by 
forced  Rayleigh  scattering  (FRS)  using  methyl  red  as  a  tracer. 

The  least-squares  fits  for  the  free- volume  approach  were  done  using  Delta 
Graph  Professional  (DeltaPoint  Inc.,  Monterey,  CA)  which  uses  the  Levenberg- 
Marquardt  method  to  fit  the  data.  For  the  kinetic  theory  approach,  Mathematica 
(Wolfram  Research,  Champaign,  IL)  was  used  for  data  fitting. 
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Results 

Diffusion  experiments  were  performed  on  the  polymer-solvent  systems  at 
25°C  and  polymer  weight  fractions  between  0  and  0.50.  The  results  of  these 
measurements  are  given  in  Table  I.  Figure  1,  shows  the  reduced  diffusion 
coefficients  of  these  systems  plotted  as  a  function  of  polymer  volume  fraction. 
Also  shown  in  Figure  1  are  the  predictions  of  the  Maxwell-Fricke  model  and  the 
Mackie-Meares  model  for  the  dependence  of  the  diffusion  coefficient  on  polymer 
volume  fraction.  We  remind  the  reader  that  plotting  the  log  (D/Dq)  has  the  effect 
of  emphasizing  the  differences  between  the  different  systems,  whereas,  plotting 
D/Do  highlights  the  similarities  in  the  different  systems.®'^® 

In  order  to  compare  the  predictive  ability  of  the  models,  best  fits  to  the  data 
for  the  free-volume  and  kinetic  theory  models  are  shown  in  Figures  2-7  along 
with  the  experimental  data.  Figure  2  shows  the  data  for  toluene,  ethylbenzene, 
and  cumene  in  PS.  These  three  solvents  are  very  similar  chemically.  The  solvent 
self-diffusion  coefficient  decreases  with  increasing  molecular  size  at  any  given 
concentration.  The  data  for  toluene  is  somewhat  lower  than  previously 
reported. A  similar  trend  is  seen  in  Figure  3  for  MMA  and  MEK  in  PMMA. 

The  diffusion  results  for  MEK,  chloroform,  and  t-butyl  acetate  in  PS  are 
shown  in  Figure  4.  Since  both  the  solvent  size  and  structure  are  varied  for  this 
group,  results  different  from  Figures  2  and  3  are  seen.  Chloroform  is  the 
smallest  molecule  of  the  three  solvents  (but  the  heaviest),  its  diffusion  coefficient 
is  much  lower  than  MEK’s  at  all  concentrations.  In  addition,  as  the  polymer 
concentration  increases  the  chloroform  data  shows  a  markedly  different 
concentration  dependence.  These  effects  are  probably  the  result  of  chloroform’s 
different  structure  and  polarity  compared  to  the  other  solvents. 
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The  effect  of  varying  the  polymer  for  a  given  solvent,  in  this  case  MEK  in 
PS,  PMMA,  or  PBMA  is  shown  in  Figure  5.  The  log  D  values  in  this  plot  have 
been  adjusted  to  spread  the  data  out  for  ease  of  viewing.  An  explanation  of  the 
adjustment  is  given  in  the  figure  caption.  This  plot  shows  that  the  pol3nner  has  a 
significant  effect  on  the  concentration  dependence  of  solvent  diffusion.  Figures  6 
and  7  contain  the  FRS  data  showing  that  both  the  kinetic  and  free-volume  theories 
are  not  limited  in  their  application  to  data  produced  by  NMR  measurements. 
Figure  6  also  shows  a  dependence  on  the  polymer  similar  to  that  observed  in 
Figure  5. 
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Discussion 

Maxwell-Fricke  modeL^^*^^  The  Maxwell-Fricke  model  is  based  on  electric 
conductivity  and  permeability.  It  assumes  immobile  impenetrable  species 
suspended  in  a  mobile  solvent  continuum.  The  Maxwell-Fricke  model  predicts 
the  maximum  reasonable  diffusion  coefficient  in  the  presence  of  obstructions 
when  the  obstructions  are  spheres. Assuming  this  theory  holds  even  for 
relatively  small  pol3aner  segments  which  are  immobile  relative  to  the  solvent,  the 
Maxwell-Fricke  equation  is 


D(l-<t»  1  -  (>• 


Do 


1^^ 

X 


(2) 


where  D  is  the  diffusion  coefficient,  Do  is  the  diffusion  coefficient  of  the  pure 
solvent,  <t)  is  the  polymer  volume  fraction,  (|)‘  is  the  volume  fraction  of  polymer  plus 
any  "non-diffusing"  solvent,  and  x  is  the  shape  factor  (ranging  from  1.5  for  rods 
to  2.0  for  spheres).  If  we  assume  <()=  <{)'  Equation  2  becomes^® 


D_ 

Do 


(3). 


Equation  3  is  the  form  of  the  Maxwell-Fricke  model  plotted  in  Figure  1.  As  can  be 
seen,  this  model  severely  overestimates  the  diffusion  coefficient  and  the  fit  is 
worse  at  higher  concentrations.  However,  if  one  allows  for  some  "non-diffusing" 
or  bound  solvent.  Equation  2  has  been  shown  to  effectively  fit  diffusion  data  as  well 
as  estimate  the  amount  of  non-diffusing  solvent.^®  This  approach  differs  from  the 
other  approaches  presented  in  this  work  in  that  in  order  to  fit  the  data  the  solvent 
needs  to  be  partitioned  into  two  fractions. 
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Mackie-Meares  model.  The  Mackie-Meares  model  is  based  on 

tortuosity  (the  increased  path  length  between  two  points  due  to  obstructions)  and 
simple  cubic  lattice  model  with  similar  size  sites  for  the  solvent  and  a  polymer 
segments. 


D_ 

Do 


(1+0)^ 


(4) 


The  Mackie-Meares  model  does  a  reasonable  job  of  fitting  the  experimental  data 
in  Figure  1.  While  to  a  first  approximation  the  data  are  fairly  well  fit,  there  are 
subtle  differences  between  the  various  systems  that  this  model  does  not  account 
for.  Thus,  even  though  this  equation  contains  a  simple  dependence  on  polymer 
concentration,  further  models  which  allow  for  system-dependent  variations  are 
required  for  accurate  modelling. 

The  Modified  Enskog  model.^®’^®  As  stated  previously  this  is  based  on 
kinetic  theory  and  takes  into  account  momentum  transfer  between  rigid  spheres. 
It  was  developed  to  describe  diffusion  of  liquids  and  gases  in  membranes,  but  the 
equations  fit  solution  diffusion  data  quite  well.  The  diffusion  coefficient  of  a 
tracer  molecule  in  a  pol3mier  solution  is  described  by  the  equation^® 


where  Dt  is  the  self-diffusion  coefficient  of  the  tracer  molecule,  DtM  is  the  mutual 
diffusion  coefficient  of  the  tracer  with  respect  to  the  polymer,  Dts  is  the  mutual 
diffusion  coefficient  of  the  tracer  with  respect  to  the  solvent,  DsM  is  the  mutual 
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diffusion  coefficient  of  the  solvent  with  respect  to  the  polymer  and  Xt  is  the  mole 
fraction  of  tracer  in  solution.  For  a  system  where  xt  is  very  small,  or,  where  the 
tracer  and  the  solvent  are  the  same,  Equation  5  reduces  to 


Dt  DtM  Dts  ■ 


(6) 


This  situation  will  hold  true  for  all  of  our  experimental  data.  For  the  FRS  data, 
the  tracer  mole  fraction  was  on  the  order  of  10'^  to  For  NMR  diffusion 

measurements,  the  second  condition  holds  where  DsM  =  Dim  The  respective 
mutual  diffusion  coefficients  are  expressed  as^® 


kT(ma  +  mp)l2  i 
27nnamp  J  ga|3(Cfa(5) 


(7) 


Here  the  subscripts  a  and  P  are  generic  subscripts  for  t,  s,  and  M  in  Equations  5 
and  6,  n  is  the  number  of  molecules  per  milliliter.  Gap  is  the  hard  sphere  collision 
radius  of  components  a  and  P  which  is  the  sum  of  the  radii  of  the  two 
components,  k  is  the  Boltzmann  constant,  T  is  the  absolute  temperature,  mj  is  the 
molecular  mass  of  component  i,  and  gap(cfap)  is  the  rigid  sphere  contact  radial 
distribution  function.  This  distribution  function  is  obtained  from  scaled  particle 
theory^® 

a  .  1  ,  6^2  [  RgRp  ]  ^  12^2  j  RgRp  ?  (8) 

1-^3  (l  -  ^3n^« 

where 

5,=|inj(2Rj)‘. 

j=l 


(9) 
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■u  is  the  number  of  components  solution.  The  j's  are  the  components  a,  (3,  or  y. 
Therefore,  nj  is  the  number  of  molecules  per  milliliter  of  component  j  and  Rj  is 
the  molecular  radius  of  component  j. 

The  diffusion  coefficient  of  component  a  with  respect  to  the  polymer 
molecule  is  DaM-  However,  the  solvent  and  tracer  molecules  do  not  interact  with 
an  entire  polymer  chain  at  one  time.  Consequently,  it  is  appropriate  to  scale  the 
diffusion  coefficient  with  respect  to  part  of  the  chain,  for  example  an  individual 
monomer  unit.  Dam 

d.m  =  again  ,10) 

where  nm  is  the  number  of  monomer  units  per  milliliter.  Substitution  of 
Equations  7  through  10  into  either  Equation  5  or  6  gives  a  rather  complicated 
dependence  of  the  tracer  diffusion  coefficient  on  concentration.  The  resulting 
equation  is  a  function  of  known  constants  (the  molecular  mass  of  each 
component,  mj),  concentration  terms  (the  number  of  molecules  per  millimeter,  n; 
the  number  of  molecules  per  millimeter  of  tracer,  ni;  the  number  of  molecules 
per  millimeter  of  solvent,  ng;  and  the  number  of  monomer  units  per 
millimeter,  nm),  and  the  molecular  radii  of  each  component  which  are  used  as 
fitted  parameters.  To  obtain  the  fits  to  the  data  shown  in  Figures  2  through  7,  a 
non-linear  least  squares  fit  to  Equation  6  was  performed  yielding  values  of  the 
radii. 

As  can  be  seen  from  the  plots,  the  Enskog  theory  does  an  excellent  job  of 
fitting  the  experimental  data  and  accounting  for  the  differences  in  solvents  which 
the  Mackie-Meares  model  did  not.  The  plots  show  that  the  modified  Enskog 
model  fits  the  data  best  at  the  lower  polymer  concentrations,  however,  at  higher 
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concentrations  in  several  systems,  the  fit  diverges  from  the  experimental  data. 
In  the  toluene/PS,  t-butyl  acetate/PS,  MEK/PMMA,  MEK/PBMA,  and  THF/PS 
systems,  the  modified  Enskog  fit  deviates  significantly  from  the  experimental 
data  at  higher  polymer  concentrations.  This  is  expected  because  it  was  assumed 
the  solvent  molecules  and  the  polymer  segments  are  rigid  spheres  that  collide. 
At  high  '^-'  polymer  concentrations,  the  polymer  chains  are  closer  and  the  solvent 
molecules  may  interact  with  more  than  one  segment  at  a  time.  In  addition,  the 
poljmier  segments  may  not  be  able  to  move  independently. 

The  solvent  diameters  and  monomer  segment  diameters  obtained  from 
these  fits  are  listed  in  Table  II.  These  results  are  all  reasonable  sizes  for  the 
solvent  molecules  and  monomer  segments,  however,  the  variation  in  size  for  PS, 
PMMA,  and  THF  for  different  systems  is  higher  than  expected.  In  the  cases  of 
PS  and  PMMA  these  variations  could  indicate  that  solvent  molecules  are 
interacting  with  more  than  one  monomer  unit  at  a  time. 

The  modified  Enskog  theory  does  not  yield  a  simple  dependence  of  diffusion 
coefficient  on  concentration.  It  does  provide  an  experimental  way  to  estimate 
solvent  molecular  size  or,  alternatively,  if  this  information  is  available  elsewhere, 
this  theory  could  be  used  as  a  predictive  method  for  diffusion  concentration 
dependence.  Table  III  contains  molecular  volumes  determined  by  group 
contribution  methods  which  could  be  used  as  a  source  of  molecular  diameters. 
Diameters  obtained  from  these  volumes  by  assuming  the  molecules  are 
spherical,  are  significantly  larger  than  the  diameters  determined  from  the  fits, 
given  in  Table  II.  None  of  the  molecules  in  question  are  nearly  spherical  in 
shape  and  this  approximation  apparently  overestimates  the  collision  contact 
radius  for  these  molecules.  Zhu  and  Macdonald®®  have  successfully  treated  the 
diffusion  of  nonspherical  solvent  molecules  in  polymer  systems. 
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Free-Volume  Theory. Some  of  the  problems  noted  with  the  kinetic 
theory  approach  —  such  as  isolated  free-volume  and  interaction  with  more  than 
one  monomer  unit  are  taken  into  account  in  a  free-volume  approach.  The 
Vrentas-Duda  free-volume  theory^®’^^  gives  the  concentration  dependence  of  the 
diffusion  coefficient  as 


D  =  Doi  ex 


-E 

exp 

0)1  Vi  +  0)2^V2 

RT 

(^^U|K2,.Tgi+T)  + 

fKi2) 

1  Y  1 

jcDj 

!(K22-Tg2+T) 

fll) 


where  Dqi  is  the  preexponential  factor,  E  is  the  energy  of  activation  for  a  solvent 
jump,  R  is  the  gas  constant,  T  is  the  absolute  temperature,  o)i  is  the  weight 
fraction  of  component  i  (1  is  for  the  solvent  and  2  is  for  the  polymer),  Vj  is  the 
specific  volume  needed  for  one  jumping  unit  of  component  i  to  make  a  jump,  ^  is 
the  ratio  of  the  volume  of  the  solvent  jumping  unit  to  the  polymer  jumping  unit 
which  is  not  necessarily  the  same  as  the  volume  of  one  monomer  segment,  Kn/y 
and  K2i-Tgi  are  solvent  free-volume  parameters,  and  K21/Y  and  K22-Tg2  are 
pol3nner  free-volume  parameters. 

For  a  pure  solvent  Equation  11  reduces  to 


Do  =  Doi  exd 


-E 

RT 


(12) 


Combining  Equations  1 1  and  12  gives 
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Assuming  that 


»(^ 


VT), 


which  should  be  true  at 


low  polymer  concentrations,  Equation  13  can  be  rearranged  to  give 


* 


Since  there  is  much  more  free-volume  per  unit  volume  in  a  solvent  than  in  a 
polymer, the  assumption  made  to  get  Equation  14  should  hold  up  to 
approximately  50%  polymer.  By  this  point  the  total  free-volume  provided  by  the 
solvent  will  have  been  reduced  to  the  point  that  the  polymer  free-volume  must  be 
considered.  The  exponent  in  Equation  14  is  the  free-volume  needed  for  a  polymer 
segment  to  jump  divided  by  free-volume  available  in  the  solvent. 

At  lower  polymer  concentrations  (<  50%)  it  is  useful  to  rewrite  Equation  14 


as 


AcOg 

(1-0)2) 


(15) 


where 


This  is  a  different  functional  form  than  those  previously  proposed,  with  A  being 
the  only  system  dependent  parameter.  Values  for  Kn/y  and  K2i-Tgi  for  many 
solvents  are  available  from  the  literature.  In  cases  where  these  parameters  are 
not  available  from  the  literature  they  can  be  calculated  from  temperature 
dependent  viscosity  data  for  the  solvent.  Such  data  is  readily  available  for  many 
systems. V2  is  the  volume  needed  per  gram  of  polymer,  for  a  polymer  jumping 
unit  to  jump  and  this  is  equivalent  to  the  polymer  specific  volume  at  0°K  which 
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can  be  calculated  from  group  contribution  methods.^  With  the  values  for  these 
parameters,  when  A  is  determined  by  fitting  Equation  15  to  experimental  data,  ^ 
can  be  determined.  Molecular  volumes  can  be  calculated  from  Vi  and  V2  ,  then 
using  the  experimentally  determined  4  the  number  of  monomer  segments  in  a 
pol5nner  jumping  unit  can  be  determined. 

In  contrast,  past  applications  of  free-volume  theory  to  the  concentration 
dependence  of  the  reduced  diffusion  coefficient^^*^®’®^’^®’^®  have  started  with  the 
free-volume  theory  of  Fujita®^  and  Doolittle.^®  This  approach  has  been  shown  to 
yield  a  linear  relationship  between  log(D/Do)  and  or  l/log(D/Do)  and 

In  some  cases,  pure  polymer  rather  than  pure  solvent  was  used  as  the 
reference  state,®^’®®  however,  in  these  instances  there  was  a  significant  difference 
between  experimental  data  and  theory  at  high  solvent  concentrations.  This  was 
attributed  to  deviations  from  free-volume  theory  at  high  solvent 
concentrations.®^’®®  While  this  approach  can  lead  to  values  for  free-volume 
parameters  for  the  system  it  does  not  give  information  on  the  size  of  the  polymer 
jumping  unit  which  we  can  get  from  the  Vrentas-Duda  theory. 

The  fits  of  Equation  15  to  the  experimental  data  are  shown  as  the  solid 
curves  in  Figures  2-7.  As  stated  above  log  D  was  plotted  rather  than  log(D/Do)  for 
ease  of  viewing.  Equation  15  appears  to  describe  the  experimental  data  well,  in 
fact  in  several  cases  better  than  the  modified  Enskog  theory.  However,  at  lower 
polymer  concentrations  this  free-volume  approach  does  not  fit  the  data  as  well  as 
the  kinetic  theory  approach.  This  because  at  very  low  polymer  concentrations 
excluded  volume  does  not  play  as  significant  a  role  as  collisions  between  the 
solvent  molecules  and  the  polymer  chains.®^’®®  Table  III  shows  the  A’s  and  4’s 
determined  from  the  fits  as  well  as  the  number  of  monomer  units  per  pol3mier 
jumping  unit. 


Vfaggoner,  Blum,  and  MacElroy  17 


The  free-volume  approach  does  give  a  relatively  simple  dependence  of 
diffusion  coefficient  on  concentration  and  accounts  for  differences  between 
systems.  Although  this  method  is  unable  to  be  completely  predictive,  the 
experimental  results  indicate  that  A  will  be  close  to  -1.  Most  of  the  experimental 
values  for  A  shown  in  Table  III  range  from  -1.2  to  -0.4.  Table  Ill  also  gives  the 
number  of  monomer  units  per  polymer  jumping  unit  determined  from  the 
experimental  fits.  For  the  most  part,  the  values  seem  reasonable,  indicating  that 
between  one  and  two  monomer  units  is  involved  in  a  jumping  unit  and  that  this 
varies  some  what  for  a  given  polymer  depending  on  the  solvent.  The  variation  of 
the  polymer  jumping  unit  for  solvents  in  the  same  polymer  might  be  consistent 
with  chain  motions  occurring  with  a  broad  range  of  torsion  angle  jumps  rather 
than  the  same  set  jump  each  time  (say  120°).  This  type  of  chain  motion  has  also 
been  observed  in  the  bulk  atactic  poly(propylene)^®  and  polystyrene^®  near  the 
glass  transition  temperature  by  deuterium  2D  exchange  NMR. 

Some  of  the  systems  studied,  yielded  jumping  units  of  less  than  one 
monomer  unit  such  as  PMMA  in  MA  (0.70),  PS  in  THF  (0.33),  and  PIP  in  THF 
(0.699).  For  longer  mers  such  as  those  in  polyisoprene  (PIP),  not  all  of  the 
backbone  atoms  need  to  be  involved  in  segmental  motions.  This  could  account  for 
the  lower  value  found  for  it.  For  PMMA,  side  chain  motions,  may  account  for  the 
jumping  unit,  possibly  explaining  the  PMMA  in  MM  A  jumping  unit  of  0.70. 
However,  the  jumping  unit  for  PMMA  in  MEK  is  1.62  monomer  units.  This 
difference  in  the  size  of  the  jumping  unit  for  the  same  polymer  may  be  due  to  the 
way  the  different  solvents  interact  with  the  polymer.  However,  these 
explanations  for  jumping  units  of  less  than  one  monomer  unit,  will  not  work  for 
PS,  with  only  two  atoms  in  the  backbone  and  an  aromatic  group  for  a  side  chain. 
In  determining  the  monomer  units  per  jumping  unit  from  it  was  assumed  that 
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an  entire  solvent  molecule  was  the  solvent  jumping  unit.  If  this  is  not  true,  then 
the  results  will  be  distorted.  This  seems  an  improbable  cause  for  the  extremely 
small  polymer  jiunping  unit  in  the  THF/PS  system.  The  data  for  this  system 
came  from  FRS  measurements  so  the  diffusion  coefficient  is  actually  that  of 
methyl  red  not  THF.  However,  since  the  methyl  red  is  at  such  a  low 
concentration  it  should  make  little  difference  in  the  free- volume  available,  it 
would  seem  that  it  would  follow  the  same  concentration  dependence  as  the 
solvent^®'^^  (which  is  responsible  for  nearly  all  the  free-volume).  However,  if  the 
methyl  red  would  preferentially  associate  with  the  polymer  it  would  show  a 
greater  concentration  dependence  and  thus  falsely  indicate  a  smaller  monomer 
imit  per  polymer  jump  unit. 

The  assumption  that  preferential  association  will  result  in  increased 
concentration  dependence  is  borne  out  in  Figure  7.  Lodge  et  al.^®  *^^  have  shown 
that  methyl  red  preferentially  binds  to  polyCvinyl  acetate)  in  a  toluene/poly(vinyl 
acetate)  solution  via  hydrogen-bonding.  The  data  shown  in  Figure  7  are, 
experimentally  determined  data  (squares)  and  the  data  with  Lodge's  correction  to 
remove  the  hydrogen-bonding  effect  (diamonds),  for  the  toluene/poly( vinyl  acetate) 
system.^®  Kinetic  auid  free-volume  theory  fits  are  shown  for  both  sets  of  data  and 
it  can  be  seen  that  neither  theory  fits  the  uncorrected  data  very  well.  Table  II 
gives  the  solvent  and  monomer  radii  determined  from  kinetic  theory  for  both  the 
corrected  and  the  uncorrected  data.  The  uncorrected  data  gives  a  substantially 
smaller  solvent  radius  and  larger  monomer  radius.  It  would  be  difficult  to  use 
this  as  a  test  for  preferential  binding.  The  free-volume  fit  results  given  in  Table 
III  are  more  helpful,  here  the  non-corrected  data  gives  a  polymer  jumping  unit  of 
0.757  monomers,  while  the  data  corrected  to  remove  hydrogen-bonding  effects 
yields  1.75  monomer  units  per  jumping  unit.  The  later  value  seems  much  more 
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reasonable  and  supports  the  conclusions  given  above  regarding  the  effect  of 
preferential  binding  of  the  tracer  on  the  free-volume  fit  results.  Therefore 
preferential  binding  of  the  methyl  red  tracer  is  likely  in  the  case  of  THF/PS  and 
possibly  in  THF/PIP  as  well. 

Diffusion  of  MEK  in  PS,  PMMA,  and  PBMA  (Figure  5),  shows  that  the 
polymer  has  a  significant  effect  on  the  concentration  dependence  of  solvent 
diffusion.  The  free-volume  results  indicate  that  as  the  polymer  jumping  unit  per 
monomer  unit  increases,  the  decrease  in  diffusion  coefficient  with  increasing 
polymer  concentration,  is  reduced.  This  can  be  seen  since  A  decreases  from 
PBMA  to  PS  as  does  the  monomer  units  per  jumping  unit.  However,  the  volume 
of  the  jumping  unit  does  not  follow  this  trend  and  neither  does  the  volume  of  the 
monomer  unit.  This  indicates  that,  as  longer  range  chain  motions  are  needed  for 
a  polymer  jumping  unit  to  jump,  the  probability  that  the  jumping  unit  will  jump, 
decreases  leaving  more  free-volume  available  to  the  solvent. 

Comparisons.  Table  IV  shows  a  comparison  the  kinetic  and  free-volume 
theories.  The  ratio  of  the  solvent  jumping  unit  volume  to  polymer  jumping  unit 
volume,  4.  was  determined  for  the  kinetic  theory  results  by  calculating  volumes 
from  the  diameters  given  in  Table  II.  The  absolute  value  of  the  difference 
between  this  value  and  ^  determined  from  the  free-volume  approach  (Table  III)  is 
shown  in  Table  IV.  The  wide  range  in  these  values  is  indicative  of  a  monomer 
unit  and  a  polymer  jumping  unit  not  being  equivalent.  Table  III  also  shows 
for  both  fits,  for  each  system.  The  free-volume  approach  gave  the  best  fit  in  seven 
of  the  twelve  systems. 
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Conclusions 

We  have  shown  that  a  theory  must  take  into  account  specific  system  effects 
in  order  to  describe  the  diffusion  of  a  solvent  in  a  polymer-solvent  mixture  in 
detail.  The  two  theories  that  take  specific  system  effects  into  account:  the  kinetic 
theory  approach;  and  the  free-volume  approach  can  fit  the  data  adequately.  Both 
the  kinetic  and  free-volume  theories  can  follow  trends  in  solvent  size  and 
chemical  nature,  as  well  as,  the  effects  of  different  polymers  on  the  diffusion  of 
one  solvent. 

Each  approach  has  its  advantages  and  disadvantages.  However,  both 
theories  fit  the  data  quite  well,  accounting  for  differences  in  solvents  and 
polymers,  and  each  provides  unique  information  about  the  system.  The  kinetic 
theory  comes  the  closest  to  predicting  diffusion  coefficients,  needing  only  the  radii 
of  the  solvent  and  monomer  unit.  The  free-volume  approach,  however,  gives  a 
relatively  simple  dependence  of  the  reduced  diffusion  coefficient  on  concentration. 
While  it  can  not  predict  the  diffusion  coefficient  at  a  given  concentration,  without 
fitting  to  experimental  results,  it  can  be  used  to  indicate  a  likely  range  for  a  given 
diffusion  coefficient  based  on  knowledge  of  Do. 
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Table  II. 

Restilts  of  the  Kinetic  Theory  Fit 

System 

Solvent  Diameter(mn) _ Monomer  diameter(nm) 

Toluene/PS 

0.572 

0.666 

Ethyl  Benzene/PS 

0.604 

0.662 

Cumene/PS 

0.641 

0.680 

t-Butyl  Acetate/PS 

0.624 

0.674 

Chloroform/PS 

0.506 

0.704 

MEK/PS 

0.520 

0.650 

MEK/PBMA 

0.523 

0.704 

MEK/PMMA 

0.522 

0.538 

MMA/PMMA 

0.566 

0.573 

THF/PS 

0.449 

0.670 

THF/PIP 

0.393 

0.573 

Toluene/PVA  corrected 

0.562 

0.579 

Toluene/PVA  raw  data 


0.478 


0.742 


Table  III.  Resvilts  of  Free- Volume  Fits 


System 

A 

Solvent 

Volume'’ 

(nm®) 

Monomer  imit 
Volume'’ 
(nm®) 

Monomer 
Units  Per 
jumping  unit® 

Toluene/PS 

-0.884 

0.781 

0.140 

0.147 

1.22 

Ethyl  Benzene/PS 

-0.855 

0.738 

0.167 

0.147 

1.54 

Cumene/PS 

-1.02 

1.20 

0.187 

0.147 

1.06 

t-Butyl  Acetate/PS 

-1.02 

0.647 

0.171 

0.147 

1.80 

Chloroform/PS 

-1.22 

0.507 

0.101 

0.147 

1.35 

MEK/PS 

-0.644 

0.477 

0.108 

0.147 

1.54 

MEK/PBMA 

-0.416 

0.308 

0.108 

0.201 

1.75 

MEK/PMMA 

-0.636 

0.509 

0.108 

0.131 

1.62 

MMA/PMMA 

-1.09 

1.583 

0.145 

0.131 

0.700 

THF/PS 

-1.06 

2.23 

0.108 

0.147 

0.330 

THF/PIP 

-0.736 

1.37 

0.108 

0.113 

0.699 

Toluene/PVA  corrected 

-0.746 

0.746 

0.140 

0.107 

1.75 

Toluene/PVA  raw  data 

-1.73 

1.73 

0.140 

0.107 

0.757 

and  the  niunber  of  monomer  units  per  jumping  unit  were  determined  from  A. 
^The  solvent  and  monomer  unit  volumes  were  calculates  from  gi  oup  contribution 
methods.®® 
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Table  IV.  Comparison  of  Free-Volume  and  Kinetic  Theory  Results 


System  1 

I4fv-4ktl 

x\t® 

best  fit'’ 

Toluene/PS 

0.147 

0.0162 

0.00307 

hr 

Ethyl  Benzene/PS 

0.022 

0.00607 

0.00268 

fv 

Cumene/PS 

0.363 

0.00469 

0.0158 

kt 

t-Butyl  Acetate/PS 

0.147 

0.0654 

0.00941 

fv 

Chloroform/PS 

0.136 

0.00240 

0.000296 

fv 

MEK/PS 

0.035 

0.000918 

0.00477 

kt 

MEK/PBMA 

0.102 

0.0139 

0.00144 

fv 

MEK/PMMA 

0.404 

0.0202 

0.00139 

fv 

MMA/PMMA 

0.619 

0.00956 

0.0428 

kt 

THF/PS 

1.927 

0.0853 

0.00998 

fv 

THF/PIP 

1.043 

0.00130 

0.00358 

kt 

Toluene/PVA  corrected 

0.168 

0.00545 

0.0248 

kt 

“The  %^'s  are  the  sinns  of  the  squares  of  the  residuals  for  the  kinetic  theoryCkt) 

and  the  free-volume  theory(fv)  for  each  fit.  '’The  last  column  indicates  which 
method  fit  the  data  best. 
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Figure  1.  Solvent  self-diffusion  coefficient  data  and  comparison  with  the 
predictions  of  the  Maxwell-Fricke  and  Mackie-Meares  models.  (♦)  toluene-PS  , 
(•)  ethylbenzene-PS,  (O)  cumene-PS,  (A)  t-butyl  acetate-PS,(n)  methyl  ethyl 
ketone-PS,  (■  )  Chloroform-PS,  and  (O)  methyl  methacrylate-PMMA. 

Figure  2.  Log  plot  of  solvent  diffusion  data  along  with  best  fits  to  free-volume  (— ) 
and  kinetic  theory  (-)  models. 

Figure  3.  Fits  to  data  for  methyl  methacrylate  and  methyl  ethyl  ketone  in 
poly(methyl  methacrylate).  Solid  lines  represent  free-volume  fits,  dashed  lines 
represent  kinetic  theory  fits. 

Figure  4.  Fits  to  experimental  data  for  non-aromatic  solvents  in  polystyrene. 
Solid  lines  represent  free-volume  fits,  dashed  lines  represent  kinetic  theory  fits. 

Figure  5.  Fits  for  MEK  in  PS,  PMMA,  and  PBMA.  The  log  D  values  have  been 
adjusted  to  spread  the  data  out  for  ease  of  viewing.  In  each  case  the  closest  tic 
mark  to  log  Dq  (i.e.  0.9  for  MEK/PBMA,  0.7  for  MEK/PMMA,  and  0.5  for  MEK/PS) 
corresponds  to  a  log  D  value  of  -8.5.  Solid  lines  represent  free-volume  fits,  dashed 
lines  represent  kinetic  theory  fits. 

Figure  6.  Fits  to  Forced  Rayleigh  Scattering  data.^®  Solid  lines  represent  free- 
volume  fits,  dashed  lines  represent  kinetic  theory  fits. 

Figure  7.  Fits  to  Lodge’s^®’^^  toluene-poly( vinyl  acetate)  FRS  data;(A)  data 
displaying-hydrogen  bonding  effect,  (■  )  data  with  the  hydrogen-bonding  effect 
removed.  Solid  lines  represent  free-volume  fits,  dashed  lines  represent  kinetic 
theory  fits. 
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predictions  of  the  Maxwell-Fricke  and  Mackie-Meares  models.  (♦)  toluene- PS  , 
(•)  ethylbenzene-PS,  (O)  cumene-PS,  (A)  t-butyl  acetate-PS,(n)  methyl  ethyl 
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Figure  2.  Log  plot  of  solvent  diffusion  data  along  with  best  fits  to  free-volume  (-) 
and  kinetic  theory  (-)  models. 
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Figure  3.  Fits  to  data  for  methyl  methacrylate  and  methyl  ethyl  ketone  in 
pol3Kmethyl  methacrylate).  Solid  lines  represent  free-volume  fits,  dashed  lines 
represent  kinetic  theory  fits. 
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Waggoner  et  al.  Figure  4. 


Figure  4,  Fits  to  experimental  data  for  non-aromatic  solvents  in  polyst3n'ene. 
Solid  lines  represent  free-volume  fits,  dashed  lines  represent  kinetic  theory  fits. 
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Figure  5.  Fits  for  MEK  in  PS,  PMMA,  and  PBMA.  The  log  D  values  have  been 
adjusted  to  spread  the  data  out  for  ease  of  viewing.  In  each  case  the  closest  tic 
mark  to  log  Do  (i.e.  0.9  for  MEK/PBMA,  0.7  for  MEK/PMMA,  and  0.5  for  MEK/PS) 
corresponds  to  a  log  D  value  of  -8.5.  Solid  lines  represent  free-volume  fits,  dashed 
lines  represent  kinetic  theory  fits. 
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Figure  6.  Fits  to  Forced  Rayleigh  Scattering  data.^®  Solid  lines  represent  free- 
volume  fits,  dashed  lines  represent  kinetic  theory  fits. 
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Figure  7.  Pits  to  Lodge’s^®’"*^  toluene-poly( vinyl  acetate)  FRS  data;(4)  data 
displaying-hydrogen  bonding  effect,  (■  )  data  with  the  hydrogen-bonding  effect 
removed.  Solid  lines  represent  free-volume  fits,  dashed  lines  represent  kinetic 
theory  fits. 


